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Abstract 

A heterogeneous brittle material characterized by a random field of local toughness 
K c {x) can be represented by an equivalent homogeneous medium of toughness, K e $. 
Homogenization refers to a process of estimating K e g from the local field K c {x). 
An approach based on a perturbative expansion of the stress intensity factor along 
a rough crack front shows the occurrence of different regimes depending on the 
correlation length of the local toughness field in the direction of crack propagation. 
A "weak pinning" regime takes place for long correlation lengths, where the effective 
toughness is the average of the local toughness. For shorter correlation lengths, a 
transition to "strong pinning" occurs leading to a much higher effective toughness, 
and characterized by a propagation regime consisting in jumps between pinning 
configurations. 
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1 Introduction 



Brittle materials such as ceramics, glass and rocks are known to be ex- 
tremely sensitive to bulk and surface defects, from which cracks can be ini- 
tiated eventually leading to failure. This extreme sensitivity calls for a sta- 
tistical analysis of crack initiation, which has been extensively develope d 
following the pioneering work of Weibull (|Weibull 19391 : Freudenthal 1968). 



This approach, describing the inception of crack propagation caused by 



initial defects dJavatilaka and Trustrum 19771 ). has been progressively ex- 



tended to account for various statistical distributions of bulk or surface 
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defect s (IMunz and Fett 19991). multiaxial criteria for critical loads on de- 
fects dBatdorf and Heinish 19781 : Evans 1 978h . inhomogeneous stress fields 



f Davies 1973j ). In this "weakest-link" approach, the analysis is focused 



on ini- 



tiation, implicitly assuming that the propagation stage is obtained systemat- 
ically over unlimited distances, due to the lower loading needed for this stage 
as compared to the initiation one. 

However, in some cases, even though cracks have been nucleated, they will 
not propagate to large distances, and their presence may still be acceptable 
in service condition of a given structure. Simple exa mples of such confined 
cracks are those induced by indentation ( Lawn 19931 ). In the latter case, the 



stress field has a rapid decay with distance from the indentation point, and 
hence a crack which can easily be nucleated may stop shortly after initiation. 
In ceramic / metal assemblies, residual stresses caused by the coefficient of 
thermal expansion mismatch can also prevent crack s to traverse the brittl e 
part so that a weakest link hypothesis does not apply ( Charles and Hild 20021 ) . 



In brittle-matrix composites, crack arrest is also observed due to the bridging 
forces induced by the fibers ( Evans 1990j ). 



Such situations require the characterization of the conditions under which a 
"macroscopic" crack may or may not propagate. The word "macroscopic" re- 
quires a specific attention. At a microscopic scale, A, a number of non-linear 
phenomena may take place in the so-called process- zone (e.g., dislocation emis- 
sion, damage initiation, or simply non-linear debonding). However, at a larger 
scale, the effect of all those confined non-linearities can be characterized by 
a toughness K c which will dictate whether a crack of stress intensity factor 
K can (when K > K c ) or cannot (when K < K c ) propagate. The situation 
considered in the present study is when, at a scale larger than A, the tough- 
ness may vary from places to places, K c (x) thus being a random field whose 
characteristics are: 

• a probability distribution function p(K c ); 

• a correlation function C defined as 

C(x) = (K c (y)K c (x + y)) y - (K c (y))l , (1) 

where the brackets {...) y denote an average over the coordinate y. 

Let us focus in the following on cases where the correlation length £ (above 
which C(x) ~ 0) can be defined, at least along the crack front. The term 
"macroscopic" refers specifically to the case of a crack whose front length 
L is much greater than £ x along the direction of the crack front. In such a 
situation, the medium can be characterized by an effective toughness that 
controls the propagation or arrest of the crack. More precisely, the local crack 
front roughness caused by the random toughness landscape can be ignored, 
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and an equivalent crack having a straight front and the same mean position 
is defined. Along this equivalent crack, a loading leading to a constant stress 
intensity factor (SIF) can be considered and the same loading will be applied 
simultaneously on the heterogeneous material. The critical SIF computed on 
the equivalent geometry which corresponds to the onset of propagation will 
thus characterize the effective toughness, K c q. 



The problem of a crack propagating in a heterogeneous toughness field has 
been considered both numerically and experimentally, for two-phase mate- 
rials having well defined geometries. I n the context of a layered toughness 
along the mean crack front, Eriksson ( Eriksson 1998j ) used an energy based 
argument to estimate the effective toughness by a rule of mixture. The by- 
passi ng of a tough inclusion has been simulated numerically in great de- 
tails ( Bower and Ortiz 1991 ; Bower and Ortiz 19931) an d corresp onding exper- 
iments have been performed ( Mower and Argon 1978j ) . Curtin ( Curtin 1997 ; 
Curtin 19 98) performed a simplified analysis of a similar problem. To ana- 
lyze the crack advance in fiber reinforced composites , instead of toughness 
distributions, strength distributions have been used (IBeverlein et al. 1997a : 
Beverlein et al. 1997bl ICurtin 19981: Landis et al. 200(3 ) . However the question 
of the large scale effective toughness was not considered. The relaxation of a 
crack di storsion due to an obstacle in a dynamic situation has also been c onsid- 
ered in ( Morrissev and Rice 2000l : Rice 200lt Woolfries and Willis 1999f ). with 
a statistical treatment addressing the question of crack front roughening. 



The heterogeneity of the toughness field induces perturbations of the 
crack front geometry. A large amount of work has been carried out 
in recent years to estimate the effect of a rough crack geometry on 
the stress intensity factor at the crack front in d ifferent cases (pla- 
nar or three-dimensiona l, static or dynamic, .. . ) dGao and Rice 19891: 
Movchan and Willis 1993 IWillis and Movchan 199^ iMovchan et al. 19981) . 
In the restricted context of a cr ack propagating in anti-plane geo metry (mode 
III), Vandembroucq and Roux ()Vandembroucq and Roux 19971 ) performed a 
second order expansion of the stress intensity factor in crack roughness. The 
second order term was shown to systematically decrease the stress intensity 
factor compared with the expected value for a straight crack. This was inter- 
preted as a "strengthening" effect, i.e., an increase of the apparent toughness 
of the material. 



Alternatively, a statistical physics approach of crack front pinning by a ran- 



tive (IPerrin and Rice 1994: 


Schmittbuhl et a 


1. 199,4 Bouchaud et al. 1993 


Ramanathan and Fisher 1997 




Ramanathan and Fisher 1998 


Tanguv et al. 19981: iBouchaud et al. 2000l Hansen and Schmittbuhl 20031) ac- 



com- 



panied by some experimental investigations (Sc hmittbuhl and Malov 1 997). 
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Most of these studies, however, focus on the statistical features of crack front 
roughness, and interesting features attached to the onset of crack propagation 
interpreted as a depinning transition. However, little attent ion has been pai d 



to the quantitative estimate of the effective toughness. In (pkoe et al. 20 02). 
the question of the statistical distribution of the macroscopic SIF at the onset 
of propagation is addressed, again revealing a universal critical behavior close 
to an effective threshold. However, the quantitative value of the latter was not 
treated. 

Crack arrest conditions have also been studied when the crack tip tra- 
verses a medium with random toughness. The problem to solve concerns 
a solid medium consisting of elastic brittle grains or potential arrest sites. 
The grain size is considered to be of constant size and toughness is con- 
stant in each grain. The toughness varies from grain to grain so that 
no spatial correlations on sc ales larger than the grain size are assumed 



(Chudnovsk v and Kunin 19871). The same hypo thesis was used by Charles et 
al. (|Charles and Hild 20021 ICharles et al. 2 003) to analyze instantaneous and 



delayed crack propagation and arrest. Jeulin ([Jeulin 19941 ) proposed a model 
in which the microstructure is assumed to be described by a Poisson mosaic. 
A Poisson tessellation defines the grain boundaries. The latter are made of 
Poisson lines in the plane for a two-dimensional medium. Instantaneous prop- 
agation and arrest conditions are investigated. 

In the following, the case of a layered local toughness, invariant along the 
direction of propagation, is considered first, since this particular texture allows 
for a closed-form answer. Then slow modulations of the toughness along the 
propagation direction are considered, in the so-called "weak pinning" regime. 
Finally, the "strong pinning" case is addressed when the correlation length of 
the toughness perpendicular to the front becomes small. In all these cases, an 
emphasis is put on an estimate of the effective toughness and comments are 
made on specific features expected in the propagation regime. 



2 Layered toughness 



A crack front parallel to the x-axis is considered, with a toughness field transla- 
tionally invariant along the propagation direction y, dK c (x,y)/dy = 0. Start- 
ing from a straight front, as the loading is increased, the crack front begins to 
develop some roughness, being still pinned in some regions and propagating 
in other parts. The front is characterized by its coordinate along the y axis: 
h(x). For simplicity, we assume that the crack remains planar, confined into 
the (x, 2/)-plane. 

The roughness of the crack front will induce a modulation of the local SIF, 
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Fig. 1. Crack propagation in a layered toughness configuration. The crack front is 
distorted such that at every point x along the front the stress intensity factor K(x) 
matches exactly the local toughness value K c (x). 

K(x), when compared to the macroscopic one K , denned with the same 
loading and a strai ght crack having the same mean position. By using a per- 
turbation analysis ( Gao and Rice 1989h . a first order solution can be derived 



K{x) =K [l 



h{x') - h(x) 



[x' — X) 



dx' 



(2) 



At the onset of propagation, the stress intensity factor has to match the local 
toughness {i.e., K(x) = K c (x)). This can be achieved through a particular 
conformation of the crack front, h*(x), which can be obtained from a mere 
Fourier transform of Eq. (2). By denoting Fourier transforms with a ~ sign, 



h*(k) 



1 K c {k) 
\k\ K. 



(3) 



The integration of Eq. (2) over x, using K(x) 
the effective toughness with K e s = K 



x 



allows one to derive 



K, 



efi 



(K c ). 



(4) 



Equation (4) is an exact result, however limited to the case of a narrow dis- 
tribution of local toughness so that the first order expansion (2) of the SIF 
modulation remains valid. We note that this simple conclusion is a general 
result when using a first order coupling term between different positions along 
the front. We note that, still using a similar linear kernel, a following section 
will show that a regime (termed "strong-pinning" ) may appear where such a 
simple conclusion breaks down. 
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-a- -fa- 
Fig. 2. A random toughness configuration is represented schematically on the left 
(a) with a grey-scale coded toughness. On the right, (b), a reference case is con- 
structed by extracting the toughness pattern along a line parallel to the propagation 
direction, and substituting to the rest a uniform toughness environment. The latter 
is determined such that the deflection Ah = h(x Q ) — (h(x)) x is zero on average over 
all accessible crack front conformations. 

3 Self-consistent homogenization 



Let us now consider the more general case of a random toughness field, vary- 
ing along the propagation direction. In contrast to the layered case, an exact 
result is not derived, but rather an approach based on a self-consistent approx- 
imation is proposed, similar to the one use d for homogenization of a randomly 



heterogeneous elastic solid (jKroener 1 967). 



A narrow strip is considered parallel to the propagation direction y, centered 
on x , and of width £ x equal to the correlation length of the toughness field in 
the x-direction. Within this strip, the random toughness is preserved. Outside 
it, one substitutes to the random toughness a homogeneous toughness, K 
as shown schematically in Figure 2. The local toughness contrast induces a 
perturbation of the equilibrium position of the crack front. Depending on the 
sign of the toughness contrast, the position h(xo) of the front at the center 
of the strip is either in front of or behind the average front position (h(x)) x . 
The value of this deflection Ah(xo) = h(xo) — (h(x)) x is at first order directly 
proportional to the toughness difference AK(xo) = K(xq)—Kq. In the process 
of homogenization, the value of K Q is chosen such that (Ah(xo)) x vanishes, 
i.e., the front shape averaged over all crack front conformations along the 
disordered strip is flat. 

Such an approach turns ou t to provide rather accurate estimates in the case 



of disordered elastic solids ((Willis 19911 ). albeit there is no way to estimate a 



priori the accuracy of the result. 

Equation (2) can be used to compute the crack front shape. The stress intensity 
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100 120 140 



Fig. 3. Schematic construction of the equilibrium position of the crack front shape: 
the intersection of the toughness profile shown as a continuous curve, with the 
coupling due to the distorsion of the crack front shown as a dotted line gives here 
a unique equilibrium position shown as a black dot. 

factor to consider is K = K all along the front except in the strip where 
K = K(x , h(xo)). Inversion of Eq. (3) leads to a logarithmic front shape, and 
thus the extension of the front along the x-direction, L, always comes into 
play in the difference Ah = h(xo) — (h). Homogeneity also dictates that the 
system size has to be scaled by the strip width, £ x . Because of linearity, Ah is 
proportional to K(x , h(x )) — K so that 



K(x , (h) + Ah) -K = 



where A is a numerical constant. 



AKn 



e. iog(£/e») 



Ah, 



(5) 



Using the linear dependence between the stress intensity factor and the dis- 
tance between the crack front within the strip and away from it, the equi- 
librium front position can be determined as shown in Fig. 3, by plotting 
K c (x , (h) + Ah) as a function of Ah together with the linear relation (see 
Eq. (5)). 



4 Weak pinning regime 



When the correlation length along the y-direction is long, so that the 
gradient dK(x,y)/dy remains always smaller than the "stiffness" S = 
AK /[£ x \og(L/£ x )], a regime referred to as "weak pinning" occurs. In this 
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case, there is a unique equilibrium position for the crack front (Fig. 3) since 
the "stiffness", S, is greater than the gradient of toughness in the propaga- 
tion (y-) direction. To compute the effective toughness K eS = K , one has to 
integrate the K c values with a measure corresponding to a uniform sampling 
of (h) . This introduces a bias in the weighting of the distribution of K c values 

K eS =\im I f Y K c {y) d -^dy 
y-^qo y Jo dy 

The first term gives the simple arithmetic average of the local toughness, sim- 
ilar to the result obtained in the layered case. The second term is proportional 
to the correlation between the local toughness and its gradient. By performing 
the change of integration variable from y to K%, it becomes apparent that in 
the case of a stationary toughness field, this second integral vanishes. Thus, 
in the weak-pinning regime, the same result as for the layered case applies 

K eS = (K c ) . (7) 

Therefore this simple result appears to be quite robust, and gives confidence 
to the self-consistent approach proposed herein. 

In this scenario, the propagation is smooth, i.e., the local crack front advance 
is a continuous function of the mean crack position. The final result concerning 
the effective toughness could als o have been der ived using different approaches, 
such as a global energy balance (Eriksson 1998), resorting to the energy release 
rate rather that the crack toughness. We however refrain from using such a 
treatment because of the following section, where this argument will be shown 
to break down. 



5 Strong pinning regime 

From the simple geometric construction leading to the equilibrium position 
of the crack front, one sees that a very different behavior appears as soon 
as the "stiffness" , S, becomes smaller than the gradient of toughness in the 
propagation (y-) direction. In this case, illustrated in Fig. 4, multiple solutions 
can be found. In generic cases, a sequence of alternatively stable and unsta- 
ble solutions is obtained. However, even ignoring those unstable positions, the 
multiplicity of solutions indicates that the position will be selected through 
the history of propagation. In the case of interest, the mean crack front po- 
sition is chosen to be monotonically increasing, so that the toughness values 
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Fig. 4. Schematic construction for the strong pinning case. The continuous curve 
shows the local toughness profile. In contrast to the previous figure, many solutions 
of equilibrium position exists as shown by the five intersections of the dotted line 
with the continuous curve (three stable positions and two unstable ones). This 
selects only a subset of accessible arrest sites shown by bold dots, separated by 
sudden jumps shown by a dashed line of slope —S. 



sampled in the strip are confined to the highest values, as if the toughness 
profile were illuminated by a grazing incidence light (i.e., exposed hill tops 
correspond to sampled equilibrium position, whereas shadowed regions corre- 
spond to unstable jumps to a new arrest position). In contrast to the weak- 
pinning case where the propagation was smooth, one observes in the strong- 
pinning case a sort of stick-slip crack propagation. This regime is the one 
addressed in the statistical ap proaches devel oped in dBouchaud et al. 1 993; 



Ramanathan and Fisher 1997 



Ramanathan and Fisher 1 



Tan guv et al. 1998: Bouchaud et al. 2000; iHansen and Schrnittbuhl 200 




To estimate K e g, it is necessary to evaluate the probability that a particular 
value of the toughness K c (x) can constitute an accessible equilibrium position. 
This implies that for all x' < x, K c (x') < K c (x) + S(x — x'). To proceed, it 
is convenient to specialize the toughness profile to a simple pattern, such as 
a piecewise constant toughness over intervals of equal size £ y (chosen equal 
to correlation length along the y direction), and without correlation from one 
interval to the next, as shown schematically in Figure (2a). 

Let P{K C ) be the cumulative toughness distribution P(K) = f p(k)dk. The 
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probability, f(K c ), that K c can be reached is written 



f(K c ) =p(K c )l[P(K c + iSQ 

i=i 

« p(AT c ) exp (l/(S£„) log(P(K c + x))dz 

This distribution is not normalized, as can be seen by the existence of shadow 
zones containing unaccessible sites. Thus, the p.d.f. of effective stable and 
accessible positions is 

and thus the effective toughness is again the simple arithmetic average of K c 
using the probability distribution <fi instead of p 

K " = WTWdk ■ (10) 

This last equation, together with the definition of / in Eq. (8), gives an es- 
timate of the effective toughness based on the self-consistent approach. Even 
though it was derived in the context of strong-pinning, in the weak pinning 
regime, Eq. (8) provides f(k) = p(k), and hence the previously derived result 
can be recovered. 

It is worth noting that the correlations along the crack front and perpendicular 
to it contribute significantly to the final formula. Consequently, the effective 
toughness is not a specific property of the material along the potential crack 
plane, but may depend on the direction of crack propagation. This result 
contrasts with energy-dissipation based arguments where such an orientation 
effect cannot appear. The fundamental basis for such a difference comes from 
the specificity of the strong pinning regime, namely the occurrence of sudden 
jumps over unstable configurations. During those jumps, the potential energy 
is transferred to kinetic energy, and this dynamic aspect is not considered 
(i.e., the implicit assumption is that most of this kinetic energy is either 
dissipated or radiated away from the crack so that it is no longer available 
for propagation). This effect explains why energy-based approaches will not 
reproduce the proposed result. Furthermore, the crack length comes into play 
in the expression of K e g in the strong-pinning regime through the particular 
dependence of the stress intensity factor on the crack front distorsion. However, 
this dependence is only logarithmic and might be difficult to observe. 

Finally, the emphasis put on the multistability of the strong pinning regime is 
com parable to parallel analyses performed in other contexts such as solid fric- 
tion ( Caroli and Nozieres 19 96). wetting phenomena (hyst eresis of wett ing an- 



gle) JaH 

amson and Gast 19981 1. or charge density waves (|Fisher 19851) where 
strong pinning is at the origin of new phenomena at the macroscopic scale as 
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compared to the microscopic one, and where equilibrium concepts and energy- 
balance arguments are no longer operational. 



6 Summary 

In this article, the homogenization of a random toughness field is addressed, 
taking into account the effect of the crack front roughness on the local stress 
intensity factor to first order in perturbation. Starting from a layered case, 
where it can be shown that the effective toughness is equal to the arithmetic 
average of the toughness, a self-consistent scheme is introduced to deal with 
more general toughness fields. This approach allows one to extend the validity 
of the effective toughness estimate to slowly varying toughness fields along 
the propagation direction, a situation referred to as weak-pinning. However, 
as the correlation length of the toughness field along the propagation direction 
decreases, a novel behavior is encountered when the crack front position is no 
longer single-valued. In the so-called strong-pinning regime, the sampling of 
the local toughness becomes inhomogeneous, and gives rise to an apparent 
strengthening. The transition between these two regimes is progressive, and is 
accompanied at the microscopic level by an unsteady propagation displaying 
sudden jumps or burst similar to stick-slip. 
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